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1 Introduction 

Diffusive shock acceleration as an astrophysical mechanism for accelerating charged 
particles has the advantage of being highly efficient. This means however that the 
theory is of necessity nonlinear: the reaction of the accelerated particles on the shock 
structure and the acceleration process must be self-consistently included in any at- 
tempt to develop a complete theory of diffusive shock acceleration. Considerable 
effort has been invested in attempting, at least partially, to do tins and it has be- 
come clear that, in general either the maximum particle energy must be restricted by 
introducing additional loss processes into the problem or the acceleration must, be 
treated as a time dependent problem (Drury. 1984). 

In stationary modified shock structures where a significant part of the down- 
stream pressure is carried by the energetic particles, because these are in part rela- 
tivistic and have a softer equation of state than a perfect monatomic gas, the com- 
pression ratio of the shock will exceed four and tend towards the limiting value of 
seven appropriate to a strong shock in a relativistic gas. If, as is almost always the 
case, the mean free path of an energetic particle and hence the diffusion coefficient 
increases with particle energy, the diffusion length scale of the particles will be much 
greater than the thickness of the shock structure above a certain energy. The acceler- 
ation of these particles can then be described by test particle theory which applied to 
a shock of compression greater than four implies a spectrum of accelerated particles 
Hatter than f(p) oc p~ 4 where p is the particle momentum. Thus without an upper 
cutoff the energetic particle pressure and energy density will diverge. 

We conclude that stationary modified shock structures can only exist for sxrong 
shocks if additional loss processes limit, the maximum energy a panicle can attain. 
This is certainly possible and if it occurs the energy loss from the shock will lead 
to much greater shock compressions. It is however equally possible that no such 
processes exist and we must then ask what sort of nonstationary shock structure 
develops. The same argument which excludes stationary structures also rules out. 
periodic .solutions and indeed any solution where the width of the shock remains 
bounded. It follows that the width of the shock must increase secularly with time 
and it is natural to examine the possibility of selfsimilar time dependent solutions. 


2 Equations 


Our basic, idea is that the upper cutoff in the particle energy spectrum is determined 
through the finite acceleration time scale by the age of the shock. From test-particle 
theory we know that the cutoff momentum p after time l is given implicitly by 
(Axford, 1981). 
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where U U+ are the upstream and downstream velocities of the plasma relative to 
the shock and « is the diffusion coefficient. Thus if the diffusion coefficient has a 
powerlaw dependence on momentum, « oc p<* with a > 0, the diffusion coefficient of 
the maximum energy particles is of order atU* and increases linearly with time^ 

If these particles carry significant amounts of energy, as strongly suggested by 
the high energy divergence of steady solutions, the associated modification of t e 
shock structure will also have a length scale which increases linearly with time. This 
suggests the possibility of selfsimilar solutions depending on the similanty vanaole 
g Jx/t. Indeed on making this ansatz and writing for the effective diffusion coef- 
ficient of the high energy particles k ( x , 1 ) = S(flt the equations ox ideal energetic 
particle hydrodynamics (see paper OG 8.1-5) reduce to the system of ordinary dif- 
ferential equations, 
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where p denotes the density, P g the pressure, U the velocity of the background plasma 
and P„ the pressure of the ultrarelativistic particles near the upper cutoff. In T e 
related problem for stationary solutions the diffusion coefficient can be eliminated by 
a change of independent variable but the adiabatic exponent of the energetic particles 
is uncertain: here k cannot be eliminated, but for the ultrarelativistic particles the 

adiabatic exponent = 4/3. , 

This system is easily found to have singularities at the points where 
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i « where the Doppler shifted velocity (note that we have transformed to an expand- 
ing coordinate system) is zero or equal to the local sound speed. A major distinction 
between such time dependent selfsimilar solutions and stationary solutions is that 
these singularities cannot be avoided; if we start with some asymptotic positive value 
of U far upstream, U - U- as £ - -oc, and integrate into the shock U Decreases 
monotonically and £ increases monotonically so that at some point V £ will f- o 
the local sound speed. Thus selfsimilar shock structures always contain a suosnock 
(whereas in the stationary case the high Mach number solutions are usually smooth,. 

The matching conditions at this subshock are essentially the same as those m 
the stationary case, namely that the energetic particle pressure be continuous and 
that the jump in the energy fiux (if any) be positive and result from injection at the 

subshock, 
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where Q is the injection energy flux. The solution downstream from the subshock 
must be spatially homogeneous, at least as far as the point where U = £ where a. 
contact discontinuity may occur. 

In the stationary case the subshock jump conditions for the plasma properties are 
simply the ordinary Rankine-Hugoniot conditions applied to a gas with ^ = 5/3 and 
the injection energy flux, though -hard to estimate, is probably sufficiently small to be 
ignored. However in the selfsimilar case the subshock is itself probably modified by 
the reaction of mildly relativistic particles and the injection flux may be significant. 

3 Results 

The above system of equations can be integrated in close analogy to the station- 
ary system. In an earlier report on this work (Beck & Drury, 1984) we considered 
the problem of constructing stationary modified structures for the subshock, but ne- 
glected injection. This was found to be possible for moderate Mach numbers, but in 
general could not be done for shocks with high Mach numbers. The problem is that 
for strong shocks the subshock, though weaker than the main shock, is still so strong 
that we return to the problem from which we started; the difficulty of constructing 
consistent modified structures for stationary shocks. The additional freedom allowed 
by including the injection flux resolves this problem, but at the expense of an ex- 
tra unknown parameter (which one can however attempt to estimate using physical 
arguments based on the concept of selfregulating injection, ef Eichler 1981, 1984). 

An interesting feature of the timedependent solutions is that the total shock 
compression is usually large, of order 10 or 20, so that both methods of resolving the 
divergence problem in strong shocks, the introduction of additional loss mechanisms 
or the explicit inclusion of secular broadening, lead to the same qualitative conclusion; 
that strong modified shocks are expected to have large compression ratios. This of 
course implies that the acceleration is efficient at accelerating particles so that, at 
least qualitatively, we have a consistent picture. 
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